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Abstract 

We give a simple tutorial introduction to the Mathematica package STRING- 
VACUA, which is designed to find vacua of string-derived or inspired four-dimensional 

= 1 supergravities. The package uses powerful algebro-geometric methods, as im- 
plemented in the free computer algebra system Singular, but requires no knowledge 
of the mathematics upon which it is based. A series of easy-to-use Mathematica 
modules are provided which can be used both in string theory and in more general 
applications requiring fast polynomial computations. The use of these modules is 
illustrated throughout with simple examples. 
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1 Introduction 

The systematic investigation of the turning points of a given potential is one of the most 
fundamental problems in physics. In particular, much of the current effort in making string 
theory reproduce the real world is concentrated in the study of effective theories with potentials 
capable of stabilising all of the moduli. The theoretical framework for deriving the potential 
given the particular scenario of compactification - be it heterotic or type 11, with or without 
fluxes - has been considerably developed. However, the analysis of the resulting effective 
theory, due to the typically complicated forms of the (super-) potential, is often computationally 
insurmountable. 

In PQ [2] a novel, algorithmic approach to this problem was presented. The key realisation of 
these papers was that the typical potential is a ratio of, albeit complicated, polynomials in the 
held variables (in the non-perturbative cases, dummy variables may be used as substitutions for 
exponential terms, recovering polynomial expressions). The most efficient method of studying 
systems of polynomials is computational algebraic geometry. 

Recent advances in computer algebra, in particular the advent of such packages as Macaulay 

2 [3] and Singular have rendered the subject of computational algebraic geometry an in- 
creasingly significant tool in string theory El E]- As outlined in [T], the starting point for 
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our analysis is a (string-derived or inspired) M = 1 supergravity theory with complex chiral su- 
perfields 0j, given in terms of the associated Kahler potential K = K{(f), 0) and superpotential 
W = W{(j)). The scalar potential V of such a theory is given by the well-known formula, 

V = expiK/M') {k^^F.F, - ^1^^!') ■ (1) 

Here Fi = diW + M~'^diKW are the F-terms, Kij = g^g^, is the Kahler metric, K^^ its inverse 
and M is the Planck mass. Our primary interest is in analysing the vacuum properties of this 
potential. 

One tool which computational algebraic geometry provides us with to help in this analysis is 
the so-called saturation expansion pLj. This method takes the complicated polynomial equations 
describing the extrema of the potential of a system, and breaks them up into a series of simpler 
systems, each describing one particular locus of turning points. Further algorithmic methods 
then make it possible to single out extrema of interest and to extract the physical properties of 
the associated vacua. Another tool, called elimination, provides constraints on the parameters 
of a model which must be satisfied for vacua of given types to exist. 

The computer package Mathematica has become an indispensable tool in modern theoreti- 
cal physics and most researchers are familiar with its syntax. The algebraic geometry system 
Singular is perhaps less known to the community interested in effective supergravity theories. 
With this in mind, we have written a package for Mathematica, called STRINGVACUA, which 
externally calls Singular where appropriate while sustaining the user in the comfortable envi- 
ronment of Mathematica. A range of Mathematica modules, from low to high level, provide both 
ease of use in string theory calculations as well as flexibility for a wider range of applications. 
In particular, the tools described in the previous paragraph are implemented in module form - 
allowing the user to access them with a few simple Mathematica commands. 

The methods of STRINGVACUA are of course far from restricted in use to problems in 
string phenomenology. Many of the modules will be of use in any problem involving polynomial 
computations. Another Mathematica package which makes use of Singular is [8]. In the ensuing 
pages we shall walk the reader through a quick and illustrative tutorial of how to use the basic 
commands of the package. The appendix contains a glossary of terms. 

Throughout this tutorial text in Serif font indicates a function or object defined in Mathemat- 
ica, while commands to be entered by the user into a Mathematica notebook are typeset as soj 

2 Getting started: a simple supergravity model 
2.1 A quick first look: SimpleMinimise 

Installation of the package is easy. Instructions, along with the download of the package itself, 
can be found at the following URL. 

http:/ /www-thphys. physics. ox. ac.uk/user/Stringvacua 

In this tutorial, we describe use of the package in a Unix environment, focussing on the 
Mathematica notebook interface although the package also works with the text interface or 
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running in emacs. The only significant differences from a Windows environment are during 
installation, for further details see the URL given above. Once you have the package installed 
start up Xll, if it is not already running, and then to load the package simply evaluate (note 
the single open quote!) 



«Stringvacua" 

A package header and links to examples and the browser help pages will be displayed, and 
the package is ready for calculations. For illustration, let us start with a very simple example 
based on one modulus T = t + ir, with Kahler potential and superpotential given by 

K = -3 log(T + T), W = a + bT + cT\ (2) 

respectively. Here, a, b and c are (real) parameters. 

We will first study the vacuum structure of this theory using the SimpleMinimise command, 
which is designed to provide an elementary illustration of the package's use and to familiarise 
the user with its output. Evaluating the command 

SimpleMinimise [] 

will result in the appearance of a dialogue box saying List of complex fields? (example: {T}). In 
our case our complex field is indeed T, thus we simply type {T} and click OK. The next two 
dialogue boxes ask for the real and imaginary parts of the complex fields. If we want these to be 
given hj T = t + ir then we enter {t} followed by {\ [Tau] } (using the Mathematica syntax). 

We are then prompted to provide the Kahler potential, -3 Log [T+Con jugate [T] ] , followed 

by the superpotential, which for our example is a + bT + cT^2 and then the parameters 
in the model, {a,b,c} . Finally (in the unix release) the package asks how much time it is 
allowed to spend on each calculation that it will perform in Singular. Let us set this to 20 
(seconds) for this example, by entering 20 . 

This data having been entered, the package then works out everything else automatically. 
An xterm will appear on your screen (by default in the upper right corner) to which Singular 
writes its protocol. The details of this protocol should not concern us for now: we will simply 
use its presence as an indication that something is "happening" in our calculation. The final 
output looks as pictured in figure [H 

Let us discuss the interpretation of this result. Clearly the first few rows of the first table, 
labelled "Model", simply repeat the data defining the model which we have input. The last 
entry in this table is the scalar potential of the system, which the package has calculated using 
the standard Af = 1 four dimensional supergravity formula ([T]). 

The second of the two tables contains the results. Each row represents an irreducible 
component of the vacuum manifold. Starting with the first row, the entry {} in the first 
column indicates that this component satisfies Ft = Re(Fr) = and Fj. = lm{FT) = (in 
notation we will use throughout) and is therefore a supersymmetric vacuum. The field values 
for this vacuum are obtained by setting the polynomials in the column entitled "0-dim. ideals" 
to zero, so for the first row they are found by solving ct^ — bt — 3a = and r = 0. 
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In[2]:= SimpleMinimise [ ] 



complex fields 


{T} 


real fields 


{t} 


imaginary fields 


{ tau} 


Kahler potential 


-3 Log [T + Con jugate [T] ] 


superpotential 


a + b T + c 


scalar potential 


2 t + 3 a (b+2 c t) +2 t (t^+tau^) +b c (s t^ + S tau^) 

6t2 



Result 



non-zero F terms 


0-dim. 


ideals 




{} 


|tau, 


- 3 a - b t H 


-ct2} 


{} 


{3 b + 


2 c t, 5 b^ 


- 4 a c + 4 tau^} 


{t} 


{tau, 


-b - c t} 





Figure 1: An example of the output of SimpleMinimise. 



Similarly, the second row corresponds to another supersymmetric vacuum. The first entry in 
the last row, {t}, indicates that this extremum of the potential has 7^ and all other F-terms 
zero (here, simply Fj. = 0) and, hence, breaks supersymmetry in the t direction. SimpleMinimise 
only provides isolated turning points of the potential in field space and so all three components 
correspond to isolated extrema (as opposed to vacua with fiat directions). The absence of a 
row starting with {t, r} shows that this model has no vacuum where supersymmetry is broken 
in the directions of both t and r. 

If no warning messages are given by the package, as is the case in this example, then 
SimpleMinimise returns all of the isolated turning points of the potential which are present for 
generic parameter values. Vacua which appear only for non-generic values of the parameters 
will be discussed in section HI 

2.2 A little more detail: sugra models in STRINGVACUA 

For more serious calculations SimpleMinimise is not the most fiexible or convenient way in which 
to perform such a decomposition of the vacuum space in STRINGVACUA. Let us repeat the 
above calculation, therefore, using some mainstream package modules. 

In this approach, we first need to define our supergravity model within Mathematica. This 
is easily done using the module CreateModel. For the above case the appropriate command 
reads 

SingleTmodel = CreateModel [{T} ,-[t} ,-[\ [Tau] }, -3 Log[T + Conjugate [T] ] J 
a + b T + c T"2, Par->{a,b,c}] ; 

In general, this module has five obligatory arguments which are, in order, the list of complex 
fields, the list of the real parts of these fields, the list of imaginary parts, the Kahler potential 
and the superpotential. In addition, we have used the option Par to tell the package about the 
model parameters a, b, c. The output of CreateModel is a model list of Mathematica assignments. 
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here called SingleTmodel, of the form keyword value. Each keyword in the model list contains 
a piece of information about the supergravity model. Examples of model keywords include 
Fields for the list of complex fields, KPot for the Kahler potential, SPot for the superpotential 
and Pot for the scalar potential. For example, the potential for the above model can be obtained 
by evaluating Pot/ . SingleTmodel which returns 

2bH + 3a{b + 2ct) + 2cH{t^ + r^) + bc{5t^ + 3r^) 

This is indeed the bosonic potential of the model ([2]), calculated using ([1]). After evaluating 
CreateModel the model list contains all of the information about the supergravity model which 
is relevant for an examination of its vacuum structure. 

Hidden from the user, the package uses techniques from algorithmic algebraic geometry in 
order to decompose the vacuum space and describe its properties [H |2] . This requires that we 
process the equations of the supergravity model into a particular form. The required operations 
are performed by the Mathematica module CalcModel. 

SingleTmodel = CalcModel [SingleTmodel] ; 

After executing CreateModel and CalcModel a model list contains over 20 pieces of data, each 
referred to by its keyword, and has essentially all of the information one would want for a 
further analysis of the model. A complete list of the physical content of a model list is given in 
Appendix \^ 

As mentioned earlier, one of the key applications of STRINGVACUA is to break up the 
equations describing the vacua of a model into smaller pieces - giving one system of equations 
for each type of vacuum. Let us see how such an analysis is performed given the model list we 
have just created. We simply evaluate 

satex = SaturationExpand [SingleTmodel] 

Once this command has been executed information will start to scroll past in the Singular 
window. While this happens a new Mathematica protocol notebook appears. Within this is 
printed the data shown in figure El The terms in the second line of this notebook appear 

Calculate 4 saturation terms : 

{1, {}, "ok"} {2, {t}, "ok"} {3, {z}, "ok"} {4, {t, z}, "ok"} 
Finished ! 

Figure 2: An example of a STRINGVACUA protocol notebook. 

successively until the calculation is finished. Typically after a few seconds, the final result is 
returned: 

{{} ^ {{Gen {r, -3a - bt + ct^}, Dim {0}, Var ^ {t,r},Par {a,b,c}}, 
{Gen ^ {36 + 2ct, 56^ - 4ac + 40^^}, Dim ^ {0}, Var ^ {t, r}. Par ^ {a, b,c}}}, (4) 
{t} {{Gen {r, -6- ct},Dim {0}, Var {t,r},Par {a,6,c}}} . 
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Each entry in this hst represents an irreducible part of the vacuum manifold of our potential, 
the package has split up the equations for an extremum of the scalar potential as desired. It is 
sometimes useful to represent results, such as the above, in table form and the package provides 
a specific command, TF, for this purpose. Evaluating 

satex//TF 

we obtain the table shown in figure [31 The results in figure [3] should be compared to those 





Dim (compl . dim . ) 


Gen (Generators) 


Par (Parameters) 


Var (Variables) 


{} 


{0} 


{i:, -3a-bt + ct^} 


{a, b, c} 


{t, z] 


{} 


{0} 


{3 b + 2 c t, 5 b^ - 4 a c + 4 i:^} 


{a, b, c} 


{t, z} 


{t} 


{0} 


{z, -b - c t} 


{a, b, c} 


{t, z} 



Figure 3: An example of the output of SaturationExpand. 

in figure [T] which were obtained by applying SimpleMinimise to the same supergravity theory. 
We see that they are the same and that the presentation is similar in the two cases. In 
figure [3] we may obtain the equations for each of the vacua by setting the polynomials in the 
column Gcn(Generators) to zero. The first column tells us which F-terms are zero in each 
vacuum. Thus, as before, we see that this model has two sets of supersymmetric vacua and one 
non-supersymmetric vacuum with Ft ^ 0. As with SimpleMinimise, if no warning messages are 
given by the package, SaturationExpand returns alloi the isolated turning points of the potential 
which appear for generic parameter values. In summary, we have found the complete isolated 
vacuum structure of the model ([2]) by performing three simple commands, namely CreateModel, 
CalcModel and SaturationExpand. 

These are, in fact, the three commands being performed during a call to SimpleMinimise. 
However, as we introduce more of the package's abilities it will be useful to have all the in- 
formation we need close to hand - this information is stored in the model lists generated by 
CalcModel and CreateModel. 

3 Filtering and finding properties of vacua 

In the cases discussed above, the equations returned by SimpleMinimise and SaturationExpand 
were so simple that they may be solved trivially to find the vacua. Any physical property 
of those vacua may subsequently be calculated using the field values obtained. It may not 
always be the case that these modules simplify the equations to such an extent. We would like, 
therefore, to have an automated method for extracting the physical properties of the vacua 
which we discover, which does not require solving the returned equations. STRINGVACUA 
provides two types of facilities towards this end, filtering and Sturm queries. 

3.1 Filtering 

As STRINGVACUA performs a computation it can discard vacua which do not meet criteria 
specified by the user. This process, which we call filtering, can greatly speed up computations. 
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It avoids costly calculations being carried out on parts of the problem which are of no physical 
interest. Many of the modules of STRING VACUA have a range of optional arguments, some 
of which deal with the task of filtering. When omitted in module calls, the filtering options 
are set to sensible default values. 

Dimension filtering 

Perhaps the simplest type of filtering offered by the package is dimension filtering. By spec- 
ifying an optional argument of the form Dim-^list of dimensions, the user can choose to filter 
out any vacua which do not have a particular dimension (the number of flat directions). In 
fact, we have already used this type of filtering in the previous section! In SimpleMinimise 
and SaturationExpand this option is set by default to Dim— >0. Thus, in the previous section, 
STRINGVACUA was only keeping vacua which are dimension zero - that is, vacua which do 
not have fiat directions. Had we wished to search for vacua with 1 or 2 fiat directions instead 
we could have evaluated the following. 

SaturationExpand [SingleTModel , Dim->{1 , 2}] 

The result, {}, tells us that there are no vacua with 1 or 2 flat directions in the model defined 
in ([2]) for generic values of the parameters. 

Field value filtering 

Another example of physically uninteresting vacua which one may want to filter out are those 
which lie in an unacceptable region of field space. For example, in a four dimensional theory 
obtained by compactifying a higher dimensional supergravity, we would wish for the size of the 
hidden dimensions to be larger than the string scale. If a particular field were to represent the 
volume of the hidden dimensions, and if our conventions were such that the string scale was set 
to 1, we may therefore only be interested in vacua where this field takes a value greater than 1. 

We can instruct modules, such as SaturationExpand, to throw away all vacua which do not 
satisfy any set of conditions, as long as those conditions can each be stated as a sign condition 
on a polynomial in the fields. Let us see how this would work, using the condition described 
above, in the following model. 

ir = -31og(r + f) -log(5 + ^), W = aS + bST + cT'^ , 
T = t + ir, S = s + ia . 

We begin by creating a model list, containing all of the relevant data of the supergravity theory, 
and add to this the ancillary mathematical structures the package requires. This is achieved, 
as in section [2], by use of the CreateModel followed by the CalcModel modules. 

modelnew=CreateModel [{T,S}, -[t,s}, {\ [Tau] , \ [Sigma] } , -3 Log[T + Conjugate [T]]^ 
|||-Log[S + Conjugate [S] ] , aS+bST+c T"2, Par->{a,b,c}] ; 
modelnew=CalcModel [modelnew] ; 

^We will encounter some of the optional arguments found in STRINGVACUA modules in the following 
sections. A comprehensive list of options for any module can be found in its Mathematica help file. 
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We now have a complete model list for the case in hand: modelnew. We now wish to apply 
SaturationExpand to this model, filtering out all vacua which do not obey the condition t > 1. 
The relevant command is 

satexmodelnew = SaturationExpand [modelnew, VarPar->{a->l , b->-l, c->l},"' 
^ RRoots->"Probab" , NRoots->l, SturmPoly->{t-l} , SturmSigns->l] ; 

This module call makes use of some new optional argument^. The argument 

VarPar ^{a l,b -l,c ^1} (6) 

sets the parameters of the model (jSD to the indicated values for the duration of this module 
call. Filtering can only be performed in cases where the parameters have been fixed; whether 
or not the required filtering conditions are met would otherwise be parameter dependent and 
would not have a unique answer. Ranges of parameter values can be scanned simply by writing 
an appropriate loop over module calls in Mathematica. The remaining options are 

RRoots "Probab", NRoots^l, SturmPoly^{t-l}, SturmSigns^l (7) 

The first two of these, in order, turn on field value filtering and only keep equations that have 
at least one real solution. The second two specify t — 1 as the polynomial whose sign is to be 
checked and state that the sign of this polynomial should be positive in any vacua to be kept in 
the output. In short, these options tell the module to only keep vacua for which the polynomial 
t — 1 takes positive values, i.e. for which t > 1 

It is worth mentioning here that the calls to Singular will, of course, take different amounts 
of time on different computers. On a modern computer the above call will finish in under 
30 seconds, however on older machines it may take longer. In the Unix release there is a 
default maximum time for Singular calculations, set at 60 seconds. This is controlled by the 
option SingularMaxTime. For example, to set the maximum time to 120 seconds for a particular 
calculation, simply add to any call the optioijfl SingularMaxTime->120 . 

The result of the above call to SaturationExpand, put in table form using satexmodelnew//TF", 
is as shown in figure |H 





"Dim ( compl . dim . ) " 


"Gen (Generators) " 


"NRoots (a 

real roots ) " 


"Par (Parameters ) " 


"Var (Variables) " 


{s, t} 


{0} 


[z, a, 

-2 + s - 2 t, 
-2 - 4 t + 3 t^} 


{2} 


{} 


{s, t, a, z) 



Figure 4: A filtered output from SaturationExpand. 

Thus only a single set of vacua exist for this model with t > 1 and they spontaneously break 
super symmetry. The first entry in the results table, {s, t} , indicates that the real parts of both 

^To find out about the options available for SaturationExpand evaluate ?SaturationExpaiid. This returns a 
brief summary of the module ending with "More...". Click on "More..." to open the Mathematica help browser 
at the relevant page. 

■^To change the SingularMaxTime setting globally, refer to ?StringvacuaDptions 
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F-terms are non-zero (while the imaginary parts vanish). The second entry tells us that this 
vacuum has no fiat directions (recall the dimension filtering is set to only Iceep isolated extrema 
by default). The third entry is the polynomials which we set to zero to obtain the equations 
describing the vacua. The fourth column tells us that these equations have 2 real solutions. As 
always, since no warning messages were given during this computation, the results in figure H] 
constitute all of the vacua of the system which satisfy our filtering constraints. 

3.2 Finding properties of vacua 

Once we have found vacua we will wish to investigate their physical properties (other than those 
determined by the filtering - which we already know!). This can be achieved in an automated 
manner using the SturmQuery module. 

As an example let us examine further the result we obtained at the end of the previous 
subsection. What else might we like to know about the vacua we have found? We know, 
from using the SaturationExpand options, that the extrema lie in a sensible part of field space. 
However, we do not yet know if they are stable - have we found saddle points, minima, or 
maxima? 

To answer this question one would usually compute the eigenvalues of the Hesse matrix 
in the vacuum and check to see if they are positive, negative or zero - but this is just a sign 
test of the type described above! Below, we will show how to formulate the relevant test by 
hand, in the interest of introducing more of the package modules. However, as stability is such 
an important issue, the package admits an option 'TestStability' for SaturationExpand, which 
automates the relevant sign tests. To illustrate, we will repeat the above SaturationExpand call, 
but adding this new option. 



SaturationExpand [modelnew, VarPar->{a->l , b->-l, c->l}, 




RRoots->"Probab" , NRoots->l 


SturmPoly->{t-l}, SturmSigns->l , 


TestStability->"Hesse"] //TF 





These options, as before, instruct the package to retain only equations which have at least 
one real solution in the allowed region of field space. The additional, final, option TestStability 
has been set to " Hesse" - this instructs the package to determine the stability of the extrema 
described by the retained equations. The output, which we have here passed straight to table 
form, is just as in figure H] (so that we recover the solutions we found previously) except that 
an additional column has been added to the table, containing 





Sturm(Sturm Query) 






{{Hessian}, {{1,-1},{2,2}},{{1,-1},{2,2}} } 





How are we to interpret this? Each entry in the column is a list (one for each set of equations 
retained). The first entry in a list, {Hesse}, simply labels the type of sign test performed - 
here it reminds us that we are testing the signs of the Hesse matrix eigenvalues. There are two 
further entries (each a list itself), because the equations have two real solutions - in general, 
then, the list would contain a label and one further list for each real solution. These entries tell 
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us about the signs of the eigenvalues at each of root of the system. In this case they are both 
identical, 

{{1,-1}, {2,2}}. (8) 

The interpretation of this is simple: at both roots, the extrema have two positive (the first 
2 goes with the +1) and two negative eigenvalues (the second 2 goes with the -1). Thus the 
non-supersymmetric vacua we have been studying are actually saddle points. 

Let us take a closer look at how the package arrived at this result. The output of the 
filtered SaturationExpand call which led to this result was called satexmodelnew. The non- 
supersymmetric vacuum we found may be extracted from this list by using its label, as given 
in the first column of figure HJ as follows, {s ,t}/ . satexmodelnew which returns 

{{Gen -> {r, ex, -2 + s - 2t, -2 - 4t + 3t^} , Dim {0}, NRoots {2}, 
Var — > {s, t, a, r}. Par {}}} 

This list contains the information in the table of figure HI in a format which the modules of 
STRINGVACUA can understand. We will now pass this information to the SturmQuery module 
to find out about the stability properties of these non-supersymmetric vacua using the sign tests 
introduced above. 

In fact, using SturmQuery we can check the signs taken by any rational functions, the signs 
taken by the determinants of any matrices of rational functions, and the signs taken by the 
eigenvalues of any matrix with real eigenvalues in the vacua we have found. In our case we 
of course wish to check the eigenvalues of the Hesse matrix. The first piece of information we 
require is the Hesse matrix itself. Fortunately this has already been calculated by the package 
and stored in the model list! The relevant keyword to access this matrix, as may be discovered 
from the help files, is Hesse. Thus evaluating Hesse/ .modelnew will return the Hesse matrix 
of the system. The result is rather lengthy and so we do not reproduce it here. 

To identify the signs of the eigenvalues of this matrix in the vacua given by ([9]) we simply 
perform the following SturmQuery call, 

SturmQuery [{s ,t}/ . satexmodelnew, TestMatEigen->Hesse/ .modelnew^ 

VarPar->{a->l, b->-l, c->l}, SturmLabel->"The Hessian", RRoots->"Probab"T| 

The first argument above is the component of the vacuum space which we wish to test, as 
detailed in ([H]). The second argument is the matrix whose eigenvalues are to be checked, as we 
just described. The third argument sets the parameters in the Hesse matrix to the correct values 
- just those that were selected in the SaturationExpand call of the previous section, from which 
([9]) was obtained. The fourth argument allows a user-defined label for the results of SturmQuery. 
Finally the fifth argument tells the module which method to use in its computations. The result 
of this call is shown in figure [5] (as earlier we have put the result in table form by appending 
//TF to the end of the result of the module call). We see that there is one additional piece of 
data in this table now, as compared to figure H) the result of the sign check we have requested. 
This is precisely the same output we found using SaturationExpand[... TestStability— "Hesse"], 
thus we recover the fact that these non-supersymmetric vacua are saddle points. A very large 
range of properties of vacua can be studied using SturmQuery. For more examples we refer the 
reader to the module's Mathematica help file. 
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Var (Variables) 


{s, t, o, z] 


Par (Parameters) 


{} 


Gen (Generators) 


[z, o, -2 + s-2t, -2-4t + 3t^} 


Dim (compl . dim . ) 


{0} 


NRoots(tt real roots) 


{2} 


Sturm (Sturm query) 


{The Hessian, {{1, -1}, {2, 2}}, {{1, -1}, {2, 2}}} 



Figure 5: A example output from SturmQuery. 

4 Constraints on parameters 

As well as finding the vacua of a given M = 1 theory, our methods may be used to identify 
constraints on the parameters of a model which are necessary for vacua of a given type to 
exist. As before we shall illustrate this with a simple example. Consider the following theory, 
obtained in [9] in the context of non-geometric compactifications of type II string theory. 



= -3 log(-z(T - f ) - log(-z(S - 5)) - 3 log(-z(f/ - f/)), 
W = ao- 3aiT + 3a2T^ - a^T^ 

+ S{-bo + 36iT - 3b2T^ + hT"") 

+ 3t/(co + (Ci + Ci - Ci)T - (C2 + C2- C2)T^ - CsT^) , 



(10) 



Here S, T and U are superfields while all other quantities are parameters. The parameters obey 
the following set of conditions, which come from consideration of tadpoles and integrability of 
Bianchi identities. We have. 



= ao&3 — 3aib2 + 3a2&i — a-^bo — 16 , 
= aoCs + ai(c2 + C2 - C2) - 02(01 + ci 



ci) - flsco , 

together with the equations below and another set given by exchanging all hats and checks. 



(11) 



= C062 - C161 + ci&i - C2&0 , = C0C2 - cl + CiCi - C2C0 , 
= Ci63 - C262 + C262 - C361 , = C3C1 - cl + C2C2 - C1C3 , 
= Co63 - C162 + C162 - C261 , = C3C0 - C2C1 + C2C1 - C1C2 , 

= C162 - C261 + C2fel - C360 , = C2C1 - C1C2 + C1C2 - C0C3 . 

As an example, we shall describe how to find further constraints on the parameters of this 
system which are necessary (and in this case sufficient) for the existence of supersymmetric 
Minkowski vacua. These vacua are obtained as solutions to the equations. 



OtW = , dsW = , duW = , W = 



(12) 



together, of course with the vanishing of the conditions on the parameters given above. This 
large and complicated set of equations contains information about what values the flux pa- 
rameters can take in order for solutions to exist, but this data is entangled with all of the 
other information about the extrema (such as field values). Our objective is to separate out 
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the information concerning the constraints on the parameters. This may be achieved using the 
package routine Elimination. 

First, we must collect the various equations describing the vacua we wish to analyse. We 
create the model list for the above system using the commands Create Mo del and CalcModel as 
in earlier sections. Let us call the resulting model list model. The equations in (fT2l) . being infor- 
mation about our supergravity model, are now stored in this model list and may be extracted, 
in a form STRINGVACUA modules will understand, by evaluating the following command. 

-ClDSPot/. model, ISPot/ .model} 

See the appendix for a description of these keywords. Next we must enter the conditions on 
the parameters given above into Mathematica. Using the usual Mathematica syntax, these are 
included as a list of the polynomials which the conditions dictate vanish, let us call it Constraints. 

Constraints={aO b3 - 3 al b2 + 3 a2 bl - a3 bO - 16, . . .}| 

Now we have entered all of the relevant equations we may simply plug them into the module 
Elimination to obtain our constraints. The first argument in this module call is simply the collec- 
tion of equations we wish to vanish. The optional argument list VarPar instructs the module to 
eliminate some set of variables from the equations it is passed. Here, we set VarPar->{T,S,U} 
in order to eliminate the field variables T, S and U from the equations. A final argument Meth 
is optional. Here we use it to select one of the faster methods of elimination available. The full 
command to be entered is 

NewConstraints = Elimination [{IDSPot/ .model , ISPot/. model, Constraints},! 
VarPar->-[T,S,U}, Meth->"Prod"] ; 

We are left with NewConstraints which contains a list of polynomials. As usual we may display 
this object in a nice form by evaluating NewConstraints//TF . The result is shown in figure 
[6l When set to zero, the polynomials in the box "Gen (Generators)" represent the constraints 
on the flux parameters which are necessary and sufficient for the existence of supersymmetric 
Minkowski vacua0. Note that in figure [6] we have relabelled the parameters in an obvious way 
for ease of input into Mathematica (ctl for ci, etc). 

Use of constraints, such as those found in this section, is one of the ways in which vacua 
which only appear for non-generic values of the parameters can be found using STRINGVACUA. 
One can simply solve for a sufficiently large number of parameters until the constraints for the 
required type of vacuum are automatically satisfied. Alternatively if one is running loops over 
calls to Saturation Expand, one for each set of fiux parameters, the constraints in this section 
can be used to greatly speed up the calculation. Any set of parameters which do not satisfy 
the constraints need not be investigated. 

''Often, one can apply one of the other STRINGVACUA modules, PrimDec, to the output of Elimination in 
order to separate this list into smaller, simpler subsets, the vanishing of each subset representing a different way 
to solve the full set of constraints. 
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Var (Variables) 


{aO, al, a2, a3, bO, bl, b2, b3, cO, c3, ccl, cc2, chl, ch2, ctl, ct2} 


Par 


Parameters) 


{} 


Gen 


Generators) 


|ct2, ctl, ch2 + ct2, chl + ctl, cc2 - ch2, 
ccl - chl, c3, cO, 16 + a3 bO - 3 a2 bl + 3 al b2 - aO b3, 
16 a3^ bO^ - 96 a2 a3 bO bl - 288 a2^ bl^ + 432 al a3 bl^ + 54 a2^ bl^ - 
81 al a2 a3 bl^ + 27 aO a3^ bl^ + 432 al a3 bO b2 - 27 a2^ a3 bO^ b2 + 
48 al a3^ bO^ b2 - 288 aO a3 bl b2 - 18 al a2 a3 bO bl b2 - 
45 aO a3^ bO bl b2 - 54 al a2^ bl^ b2 + 81 al^ a3 bl^ b2 - 

27 aO a2 a3 bl^ b2 + 54 aO a2 a3 bO b2^ + 27 aO al a3 bl b2^ - 27 aO^ a3 b2^ - 

288 al a2 bO b3 - 32 aO a3 bO b3 + 27 a2^ bO^ b3 - 45 al a2 a3 bO^ b3 + 

432 aO a2 bl b3 - 27 al a2^ bO bl b3 + 54 al^ a3 bO bl b3 + 

48 aO a2 a3 bO bl b3 + 18 aO a2^ bl^ b3 - 81 aO al a3 bl^ b3 - 

144 aO al b2 b3 + 27 al^ a2 bO b2 b3 - 54 aO a2^ bO b2 b3 - 

51 aO al a3 bO b2 b3 + 27 aO al a2 bl b2 b3 + 45 aO^ a3 bl b2 b3 - 

27 aO al^ b2^ b3 + 27 aO^ a2 b2^ b3 + 16 aO^ b3^ - 27 al^ bO b3^ + 

45 aO al a2 bO b3^ + 27 aO al^ bl b3^ - 48 aO^ a2 bl b3^ + 3 aO^ al b2 b3^} 



Figure 6: Constraints on parameters obtained using Elimination. 



5 Non-perturbative effects 

So far we have addressed models for which the scalar potential, the F-terms, and all other 
relevant quantities are all rational functions of the fields. In many cases, such as in models 
where non-perturbative contributions to the superpotential play an important role, we may be 
interested in more complicated situations. These can be addressed using the STRINGVACUA 
package by the simple expedient of replacing any non-polynomial terms in the equations with 
dummy variables. As in earlier sections we will illustrate this with some examples. 

5.1 Creating non-perturbative models 

Consider the following toy model as an example. 

W = T + exp{-S), K = -\og{T + T)-3\og{S + S) , ^^^^ 
T = t + ir , S = s + ia 

The supergravity model is set up in the package using CreateModel, just as before. 

model = CreateModel [{T,S}, {t,s}, {\ [Tau] , \ [Sigma] } , 
^ -Log[T + Conjugate [T]] - 3 Log[S + Conjugate [S] ] , T + Exp[-S]] ; 

Before we can apply our methods we must find a way to transform the system into a set 
of equations STRINGVACUA can understand - essentially into polynomials. This may be 
achieved by using dummy variables to represent non-algebraic objects. A good choice of dummy 
variables, in this case, is the set. 



X = exp(— S"), F = exp(— 5") , x = exp(— s), y = exp{—ia) 
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(14) 



These variables are built into the model as a set of replacements included in CalcModel (which, 
recall, takes supergravity objects and converts them into expressions the package can use) 
through two options. These are CDummies— a list of replacements for the superfields, and 
Dummies^..., a list of replacements for the real fields. The dummy variables defined in ([1] 
are entered as shown below. 

fullmodel = CalcModel [model , CDummies->{Exp [-S] ->X, Exp [-Conjugate [S] ] ->Y} J 
Duinmies->{Exp [-s] ->x , Exp [-I\ [Sigma] ] ->y}] ; 

To see the effect of the dummy variables, compare the supergravity potential. 

Pot/. fullmodel 



e-2^ 



3 + 12s + As' + 3e'H' + 4e^V^ + 6e"(2s - l)t cos(a) - 6e'(2a + l)r sin(a) , (15) 



A8sH 

with its numerator, converted into a polynomial, 
IPot/. fullmodel 

{Var {t, s, r, a, x, y}, Par {}, Gen {—3tx + 6stx + St'^y + Sx'^y + 12sx'^y + As'x'^y 
— 3txy^ + Gstxy'^ + 3ixT + Gisxr — Sixy'^r — Gisxy'r + Syr^}} . 

(16) 

Here, the dummy variables {x,y} have been appended to the list of true variables {t, s,r, cr} 
and the potential has been reduced to a polynomial as desired. 



5.2 Finding vacua 

The routine SaturationExpand may be used to identify and classify all vacua of this model in a 
few seconds. Note that we are working in an enlarged field space with dummy variables, in which 
true isolated vacua may appear to have fiat directions. To ensure that this is accounted for, we 
adjust the Dim option in SaturationExpand to retain vacua with any number of fiat directions 
up to and including the number of (real field) dummy variables included. SaturationExpand is 
then called using 

satexNP = SaturationExpand [fullmodel , Dim ->{0,1,2}] 

which returns a list satexNP of vacua classified by their SUSY breaking in a manner which is 
now familiar (see, for example, figure H]). Consider the vacuum space in which the imaginary 
parts of both F-terms vanish but the real parts are non-zero. We can extract the vacua of 



this kind from the full list by using their identifier, thus 


{s,t}/. satexNP 


. This list contains 


two potential vacua, as seen by evaluating 


Length [{s , t}/ . satexNP] 


. Extracting one of these 



pieces using ({s ,t}/ . satexNP) [ [1] ] we find 

{Gen ^ {r, 1 + y, -3t^ + 7tx + 2stx + 2x^ + Asx'^, -7t - 2st + x + 8sx + 4s^x, 

- Ut^ + 2st^ + lltx + + 16sx^ -15t - 2st + 2sH - 3x - 4sx, 

3t^ - 2At'x + 23tx2 + lOx^ + 2Asx^, 3t^ - 18t^x + llt^x^ - 28tx^ - 4a;^ 

- 18 - 78s - 41s2 + 4s^ + 4s^}, Dim ^ {2}, Var ^ {s, t, x, y, a, r}. Par ^ {}} 
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Thus, even in these non-perturbative cases, Saturation Expand can be used to break up the 
equations describing the vacuum space into smaller, more manageable pieces. In the above we 
could easily solve the equation given by setting the last polynomial in Gen to zero to find s and 
hence x. The solutions for r and y (and thus a) are evident from the first two polynomials. 
These solutions could then be plugged into the remaining polynomials to find t. In other words, 
in many cases the decomposition of the vacuum space provided by SaturationExpand makes the 
equations to be solved trivial. In situations where this is not the case the algorithm developed 
in ^ should be used. This algorithm can be executed using the modules of the package. This 
is, however, more involved than the level of this simple tutorial and we refer the user to [2] and 
to the package help files for more information. 

The essential point we wish to convey is that one of the core keys to the power of STRING- 
VACUA - breaking up the equations for the vacuum space into smaller pieces - is applicable 
even in these non-polynomial cases. 

5.3 Constraints in non-perturbative models 

The other application of STRINGVACUA which we have introduced in this tutorial, finding 
constraints on parameters, also works in the case where non-perturbative effects are present. 
As an example consider the supergravity model, 

K = -Slog (T + f ) - 3 log (Z + Z)- log (S + S) , 

W = i{^ + ieT) + (e + ipT)Z + i(^ + iqT)Z'^ /2 + A exp(-c5) , (18) 
T = t + ir, Z = z + i(, S = s + ia . 

which is obtained by compactifying heterotic string theory on a generalised half-flat manifold in 
the presence of gaugino condensation pHl [TT] . We can input this model into STRINGVACUA 
using CreateModel and CalcModel as usual, using the CDummies and Dummies options to specify 
the replacements for the exponentials. To clarify the presentation we here first define the Kahler 
and superpotentials. 



K = -3 Log[T+Conjugate[T]] -3 Log [Z+Conjugate tZ] ] -Log [S+Con jugate [S] ] ; 


W = I(\[Xi]+I e T) + (\[Epsilon]+I p T) Z + I (\[Mu]+I q T) Z~2/2 




+\ [Lambda] Exp [-c S] ; 







These are then entered into CreateModel using the identifiers K and W to create the supergravity 
model, 

NPmodel = CreateModel [{T,S,Z}, it,s,z}, {\ [Tau] ,\ [Sigma] ,\ [Zeta] }, K, wj 
Par->{\[Xi], e, \[Epsilon], p, \[Mu], q, \ [Lambda] , c}] ; 

Now we use CalcModel to create the extended model list, including our choice of dummy vari- 
ables for the exponential function. 



'NPmodel = CalcModel [NPmodel , CDummies->{Exp [-c S]->X, Exp[-c Conjugate [S] ] ->Y} , 


Dummies->{Exp [-C s] -> x, Exp[-1 c \ [Sigma] ]-> y}] ; 
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There is one condition on the parameters which follows from the integrability of the Bianchi 
identity, qe — pfi = 0, which we enter simply as 

^Icondition = q \[Epsilon] - p \ [Mu] ; 

We can find further constraints on the parameters of this model which are necessary for the 
existence of supersymmetric vacua satisfying certain physical conditions. For various physical 
reasons we might be interested in vacua with field values t = 1, s = 24, and a = 0, for example. 
We require, then, that the F-terms, and the polynomials t — 1, s — 24 and a all vanish in our 
vacua. 

Note that setting the value of cr = implies a fixed value for y = exp{—ica) = 1, so we 
should also set the polynomial y — 1 to vanish. At this stage, even though we have chosen a 
value for s, the value of the dummy variable x still depends on the parameter c. As such, we 
cannot fix its value. To find our constraints, we simply collect all of our equations and then 
eliminate all of the fields using the command Elimination, exactly as in section HI 

ConstraintsNP = Elimination [{IFTerms/.NPmodel, BIcondition, t-1, s-24, 
XESigma], y-1}, VarPar-> {s, t, y, z, \ [Zeta] , \[Sigma], \[Tau]}]; 

The output of this module call is shown in figure [71 put in table form by evaluating 
ConstraintsNP//TF . 



Var (Variables) {c, e, p, q, x, e. A, /J, 

Par (Parameters) {} 

Gen (Generators) {qe-p^, e q [U - , ep/U-pq<ff, IJ - e q £^ , 

-e^; + 96ce/U + q^-96cqif, e p e - -ee + 96cee+pi£;"-95cpi|", 

-^e^lJ + 3exXiJ + 4eq^-3qxX^, e;i-xA^; + 24cxA^-q^, 

-4e^e + 3exeA + 4epi£;"-3pxA.f, ee-xeA + 24cxeA-pi;", 

-qxA + 24cqxA, -pxA + 24cpxA, xAi£;"-120cxA<f + 2304c^xAi£;", 

2epe-3pxeA-12exA/U+9x^A^/j-2p^<5f+5qxAi^, 

-p2 xA-2eqxA + 3qx2A2, 4exA-96cexA-3x2A2 + 72cx2A2, 

2exA-112cexA+1536 c2exA-x2A2 + 24cx2A2} 

Figure 7: Constraints in a non-perturbative example using Elimination. 

The returned polynomials in the entry labelled "Gen (Generators)" should be set to zero 
to obtain the desired constraints. The vanishing of these polynomials describes necessary con- 
straints on the flux parameters of this model (necessary and sufficient constraints only being 
obtained in the case of perturbative supersymmetric Minkowski vacua) for the existence of 
vacua of the desired type [H [2]. As mentioned in section IH constraints such as these can fre- 
quently be simplified using other STRINGVACUA modules - for example PrimDec. The user 
is directed to the Mathematica help files for further information. 
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6 STRINGVACUA: other applications 



Since the package is tailored to the fast manipulations of systems of polynomials, its appli- 
cations clearly extend far beyond finding the vacua of string-inspired or derived models as 
described above. Problems such as finding the dimension of the space defined by a system of 
polynomials, or how many real roots there are, etc., are ubiquitous in countless fields of modern 
research. The convenient Mathematica interface we have developed, which requires neither an 
understanding of computational algebraic geometry nor of the program Singular, should prove 
useful to researchers in diverse fields of interest. In this short section we will illustrate the 
solution of questions one may encounter in more abstract situations. 

Suppose we are given a formidable-looking system of polynomials, given in figure [Hand called 
sys, and we wish to check whether its zero- locus has dimension zero (i.e., no flat directions) 
and if so find how many real roots this system has. A system of polynomials, in the language 
of algebraic geometry, is called an ideal. Ideals are stored within STRINGVACUA as a list of 
assignments, called an ideal list, as follows. 

ideal = {Gen->{polynomials}- , Var->-[variables} , Par->-[parameters}} 

Additional assignments can be added to ideal lists by some of the STRINGVACUA modules - 
we shall see an example of this shortly. 

sys = {-7 + 7 + 35 X - 2 + 2 + 10 w - 2 w y + 51 y - X y + 

67 X y + 19 w x^ y + 16 x^ y - 4 y^ + 74 y^ + 43 w X y^ + 24 x^ y^ + 44 w y^ + 40 X y^ + 
32 y^ - 8 w z + 36 z - 10 X z + 92 X z + 58 w x^ z + 12 x^ z - 4 y z + 75 y z + 
170 wxyz + 88x^yz + 110 w y^ z + 137 x y^ z + 36 y^ z - 7 z^ + 95 z^ + 123 w x z^ + 
76 x^ z^ + 89 wy z^ + 98 xy z^ + 92 y^ z^ + 92 w z^ + 122 x z^ + 39 y z^ + 56 z^ , 

30 w'' + 95 X + 134 x^ + 83 w x^ + 18 x* + 41 y + 133 x y + 153 w x^ y + 93 x^ y + 
101 y^ + 200 w X y^ + 137 x^ y^ + 52 w y^ + 132 x y^ + 70 y^ + 21 z + 96 x z + 127 w x^ z + 
71 x^ z + 98 y z + 186 w x y z + 141 x^ y z + 76 w y^ z + 209 x y^ z + 113 y^ z + 102 z^ + 
172 wx z^ + 127 x^ z^ + 95 wy z^ + 218 xy z^ + 198 y^ z^ + 36 w z^ + 99 x z^ + 126 y z^ + 81 z\ 

21 X + 21 w x^ + 3 x^ + 63 y + 93 w X y + 30 x^ y + 90 w y^ + 72 X y^ + 27 y^ + 35 z + 
59wxz + 32x^ z + 122 wyz + 140xyz + 87y^z + 40wz^ + 57xz^ + 76yz^ + 20z^, 

9 + 34 X + 79 x^ + 54 w x^ + 24 x^ + 91 y + 192 x y + 193 w x^ y + 94 x^ y + 
133 y^ + 230 w x y^ + 148 x^ y^ + 103 w y^ + 106 x y^ + 28 y'' + 14 z + 67 x z + 89 w x^ z + 
70 x^ z + 127 y z + 236 w x y z + 158 x^ y z + 107 w y^ z + 136 x y^ z + 56 y^ z + 95 z^ + 
119 wx z^ + 126 x^ z^ + 211 wy z^ + 243 x y z^ + 99 y^ z^ + 54 w z^ + 117 x z^ + 72 y z^ + 81 z'' } 

Figure 8: A given system of polynomial equations. 

Given this general structure, the dimension in our current example is easily calculated using 
the command 

Dim/ . Dimldeal [{Var-> {x, y, z, w}, Gen-> sys}P 

As described above, the Var option specifies that the variables are {x,y, z,w}, the Gen — > 
sys dictates that the system to be studied is our list sys of 4 polynomials in 4 variables and 
there is no Par assignment in this case as there are no parameters in our system. The command 
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Dimldeal[ ] then tells us the dimension of the space thus defined; it adds the assignment Dim 
— s> {the dimension} to the ideal list. Here we have extracted the dimension from this list using 
Dim /. The answer is very quickly found to be 0. 

We can now try to find the real roots. The Mathematica command Solve[ ] is readily 
seen to be ineffective in solving this system. However, methods of computational algebraic 
geometry can first reduce this system into much simpler pieces by a technique known as primary 
decomposition. What this means is that the package can split sys up into a finite number of 
simpler components, the union of whose zero-loci is the zero-locus of sys itself - in other words, 
solving any of the simpler sets of equations gives a solution to the full set sys: 

'p^ces^ PrimDec [{Var->{x,y,z,w}, Gen->sys}] ; 

The output is rather complicated so we will not present it here; suffice it to say that pieces is a 
list of length 11, check using Length [pieces] , meaning that there are 11 simple components to 
sys, each of dimension and each represented by its own ideal listEl. Now, we can find the number 
of real roots of any 0-dimensional system of polynomials with the command NumRoots[ ]. This 
adds the assignment NRoots {number of roots} to an ideal list. We can then extract the 
number with the keyword NRoots in the same way as we extracted the dimension above. 
We can perform this for each of the 11 components by entering 

Table [NRoots/. NumRoots [pieces [[i]]] , {i, 1, 11}] 

giving us 

{{1}, {0}, {0}, {1}, {0}, {0}, {1}, {0}, {0}, {0}, {0}} 

Therefore, we see that eight of the components do not have real roots while three of them have 
a single real root. In other words, sys has a total of 3 real roots. 

We can, of course, do more. Since we have the explicit expressions for the actual pieces 
we can find what the roots are explicitly. We see from the above that the 1st, 4th and 7th 
components have real roots. Executing pieces [[1]] and similarly extracting pieces[[4]] and 
pieces[[7]] will give us the relevant equations. In fact, here they all turn out to be 

{Gen {z, y, x, w}, Dim — > {0}, Var {w, x, y, z}, Par — {}} 

We see that the Gen keyword is set to {z,y,x,w}. To find the real roots we must therefore 
solve the equations obtained by setting each of the polynomials in this list to zero - which is 
just the origin. In other words, all of the roots here are at (0, 0, 0, 0) and therefore the original 
system, sys, has a triple root at the origin and no other real roots. 

7 Summary 

We have briefiy described, illustrating our discussion with simple examples, how the Mathemat- 
ica package STRINGVACUA may be used to aid in the investigation of vacuum configurations 

'''By default, PrimDec retains only dimension pieces. As before this may be adjusted using the Dim^ 
option. 
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of models of string phenomenology. We have focused on two main applications. First, we de- 
scribed how to split the large equations describing the vacuum space of a model up into smaller 
systems, each of which can be dealt with more easily. Second, we discussed how to obtain 
constraints on the parameters in a theory necessary for the existence of given types of vacua. 

In fact, in this short note, we have only used 7 of the package's 17 modules and we have 
far from exhausted the flexibility provided by optional arguments even in those. The message 
is that there is much more STRINGVACUA can do to aid the user than the simple automated 
routines described in this document. Many of the methods imported to Mathematica from 
Singular by this package are not restricted to use in string theory, but will benefit researchers 
in any area where a problem involves polynomial computations. The papers [H [2] , along with 
the Mathematica help files for the various modules, should go a long way towards helping the 
user in apply the techniques of STRINGVACUA in a more flexible manner. This, as always 
with a computer algebra system, can lead to huge improvements in what may be achieved using 
the package. 
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A Glossary — supergravity model list keywords 
A.l CreateModel 

The command model=CreateModel [ . . . ] takes a list of fields, a Kahler potential and a su- 
perpotential, and creates a model list, here called model. A model list contains essentially all 
of the information about a supergravity theory which is useful in an analysis of its vacuum 
space. Objects, such as the F-terms and the scalar potential, can be extracted via 'keywords' 
by executing a command of the form 

keyword/ . model 

where keyword can be any of the following: 
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Keyword 


Description 


Fields 


A list of the superfields in tlie model, e.g. {T}. 


ReFields 


A list of the real parts of the superfields in the model, e.g. {t}. 


ImPields 


A list of the imaginary parts of the superfields in the model, e.g. {r}. 


ComplexToReal 


A list of the relations between the Fields and their real and imaginary 




parts, in the form {T — t + it}. 


Par 


A list of the flux parameters in the model, e.g. {a}. 


KPot 


The Kahler potential. 


SPot 


The superpotential. 


DeltaPot 


The non-supersymmetric extension to the scalar potential. 


DSPot 


A list of first derivatives of the superpotential with respect to the 




superfields, in the form "{field — > derivative with respect to field}" . 


FTerms 


A list of the F-terms, in the form "{field —>■ F-term}" . 


Pi J. 

Pot 


ihe supergravity scalar potential, calculated from the superpotential. 




Kahler potential and non-SUSY extension. 


DPot 


A list of the first derivatives of the scalar potential, in the same format 




as DSPot, above. 


Hesse 


The Hessian matrix of the scalar potential, computed with respect to 




the real fields. 



A.2 CalcModel 

CalcModel creates an extended model list which contains the keywords above and added assign- 
ments where the supergravity objects have been processed into a form suitable for our methods 
- an "I" is prefixed to the above keywords to indicate this. For example. Pot/. model is now 
accompanied by I Pot/, model. All data is accessed in the manner described above. The following 
optional inputs of CalcModel are also keywords: 



Keyword 


Description 


CDummies 

Dummies 

RDummies 


A list of replacement rules for non-algebraic functions of superfields, 
in the form "function— >^ dummy variable" . 

A list of replacement rules for non-algebraic functions of real fields, 
in the form "function^ dummy variable" . 

A list of rules for replacing non-algebraic functions with dummy vari- 
ables, in the form "function— >^ dummy variable" . 
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